Candida antarctica lipase B performance in organic solvent at varying water activities studied by molecular dynamics simulations.

Applications of lipases in low-water environments are found across a broad range of industries, including the pharmaceutical and oleochemical sectors. This includes condensation reactions in organic solvents where the enzyme activity has been found to depend strongly on both the solvent and the water activity (aw). Despite several experimental and computational studies, knowledge is largely empirical, and a general predictive approach is much needed. To close this gap, we chose native Candida antarctica lipase B (CALB) and two mutants thereof and used molecular dynamics (MD) simulations to gain a molecular understanding of the effect of aw on the specific activity of CALB in hexane. Based on the simulations, we propose four criteria to understand the performance of CALB in organic media, which is supported by enzyme kinetics experiments. First, the lipase must be stable in the organic solvent, which was the case for native CALB and the two mutants studied here. Secondly, water clusters that form and grow close to the active site must not block the path of substrate molecules into the active site. Thirdly, the lipase's lid must not cover the active site. Finally, mutations and changes in aw must not disrupt the geometry of the active site. We show that mutating specific residues close to the active site can hinder water cluster formation and growth, making the lipase resistant to changes in aw. Our computational screening criteria could potentially be used to screen in-silico designed variants, so only promising candidates could be pushed forward to characterisation.

Substrates and Cyclic Peptide Inhibitors of the Oligonucleotide-Activated Sirtuin 7.

The sirtuins are NAD+ -dependent lysine deacylases, comprising seven isoforms (SIRT1-7) in humans, which are involved in the regulation of a plethora of biological processes, including gene expression and metabolism. The sirtuins share a common hydrolytic mechanism but display preferences for different ϵ-N-acyllysine substrates. SIRT7 deacetylates targets in nuclei and nucleoli but remains one of the lesser studied of the seven isoforms, in part due to a lack of chemical tools to specifically probe SIRT7 activity. Here we expressed SIRT7 and, using small-angle X-ray scattering, reveal SIRT7 to be a monomeric enzyme with a low degree of globular flexibility in solution. We developed a fluorogenic assay for investigation of the substrate preferences of SIRT7 and to evaluate compounds that modulate its activity. We report several mechanism-based SIRT7 inhibitors as well as de novo cyclic peptide inhibitors selected from mRNA-display library screening that exhibit selectivity for SIRT7 over other sirtuin isoforms, stabilize SIRT7 in cells, and cause an increase in the acetylation of H3 K18.

Combining Molecular Dynamics Simulations and Biophysical Characterization to Investigate Protein-Specific Excipient Effects on Reteplase during Freeze Drying.

We performed molecular dynamics simulations of Reteplase in the presence of different excipients to study the stabilizing mechanisms and to identify the role of excipients during freeze drying. To simulate the freeze-drying process, we divided the process into five distinct steps: (i) protein-excipient formulations at room temperature, (ii) the ice-growth process, (iii)-(iv) the partially solvated and fully dried formulations, and (v) the reconstitution. Furthermore, coarse-grained (CG) simulations were employed to explore the protein-aggregation process in the presence of arginine. By using a coarse-grained representation, we could observe the collective behavior and interactions between protein molecules during the aggregation process. The CG simulations revealed that the presence of arginine prevented intermolecular interactions of the catalytic domain of Reteplase, thus reducing the aggregation propensity. This suggests that arginine played a stabilizing role by interacting with protein-specific regions. From the freeze-drying simulations, we could identify several protein-specific events: (i) collapse of the domain structure, (ii) recovery of the drying-induced damages during reconstitution, and (iii) stabilization of the local aggregation-prone region via direct interactions with excipients. Complementary to the simulations, we employed nanoDSF, size-exclusion chromatography, and CD spectroscopy to investigate the effect of the freeze-drying process on the protein structure and stability.

Chemometrics in Protein Formulation: Stability Governed by Repulsion and Protein Unfolding.

Therapeutic proteins can be challenging to develop due to their complexity and the requirement of an acceptable formulation to ensure patient safety and efficacy. To date, there is no universal formulation development strategy that can identify optimal formulation conditions for all types of proteins in a fast and reliable manner. In this work, high-throughput characterization, employing a toolbox of five techniques, was performed on 14 structurally different proteins formulated in 6 different buffer conditions and in the presence of 4 different excipients. Multivariate data analysis and chemometrics were used to analyze the data in an unbiased way. First, observed changes in stability were primarily determined by the individual protein. Second, pH and ionic strength are the two most important factors determining the physical stability of proteins, where there exists a significant statistical interaction between protein and pH/ionic strength. Additionally, we developed prediction methods by partial least-squares regression. Colloidal stability indicators are important for prediction of real-time stability, while conformational stability indicators are important for prediction of stability under accelerated stress conditions at 40 °C. In order to predict real-time storage stability, protein-protein repulsion and the initial monomer fraction are the most important properties to monitor.

Structural characterization of human tryptophan hydroxylase 2 reveals that L-Phe is superior to L-Trp as the regulatory domain ligand.

Tryptophan hydroxylase 2 (TPH2) catalyzes the rate-limiting step in serotonin biosynthesis in the brain. Consequently, regulation of TPH2 is relevant for serotonin-related diseases, yet the regulatory mechanism of TPH2 is poorly understood and structural and dynamical insights are missing. We use NMR spectroscopy to determine the structure of a 47 N-terminally truncated variant of the regulatory domain (RD) dimer of human TPH2 in complex with L-Phe, and show that L-Phe is the superior RD ligand compared with the natural substrate, L-Trp. Using cryo-EM, we obtain a low-resolution structure of a similarly truncated variant of the complete tetrameric enzyme with dimerized RDs. The cryo-EM two-dimensional (2D) class averages additionally indicate that the RDs are dynamic in the tetramer and likely exist in a monomer-dimer equilibrium. Our results provide structural information on the RD as an isolated domain and in the TPH2 tetramer, which will facilitate future elucidation of TPH2's regulatory mechanism.

Templated Enzymatic Synthesis of δ-Cyclodextrin.

While α-, ß-, and γ-cyclodextrin (CD) are ubiquitous hosts employed by supramolecular chemists, δ-CD (formed from nine α-1,4-linked glucopyranose units) has received very little attention. α-, ß-, and γ-CD are the major products of the enzymatic breakdown of starch by cyclodextrin glucanotransferase (CGTase), but δ-CD forms only transiently in this reaction, as a minor component of a complex mixture of linear and cyclic glucans. In this work, we show how δ-CD can be synthesized in unprecedented yields by employing a bolaamphiphile template in an enzyme-mediated dynamic combinatorial library of cyclodextrins. NMR spectroscopy studies revealed that δ-CD can thread up to three bolaamphiphiles forming [2]-, [3]-, or [4]-pseudorotaxanes, depending on the size of the hydrophilic headgroup and the length of the alkyl chain axle. Threading of the first bolaamphiphile occurs in fast exchange on the NMR chemical shift time scale, while subsequent threading occurs in slow exchange. To extract quantitative information for 1:2 and 1:3 binding events occurring in mixed exchange regimes, we derived equations for nonlinear curve fitting that take into consideration both the chemical shift changes for species in fast exchange and the integrals for species in slow exchange to determine Ka1, Ka2, and Ka3. Template T1 could be used to direct the enzymatic synthesis of δ-CD due to the cooperative formation of a 1:2 complex─the [3]-pseudorotaxane δ-CD·T12. Importantly, T1 is recyclable. It can be readily recovered from the enzymatic reaction by precipitation and reused in subsequent syntheses enabling preparative-scale synthesis of δ-CD.

In-silico study of the interactions between acylated glucagon like-peptide-1 analogues and the native receptor.

We have performed a series of multiple molecular dynamics (MD) simulations of glucagon-like peptide-1 (GLP-1) and acylated GLP-1 analogues in complex with the endogenous receptor (GLP-1R) to obtain a molecular understanding of how fatty acid (FA) chain structure, acylation position on the peptide, and presence of a linker affect the binding. MD simulations were analysed to extract heatmaps of receptor-peptide interaction patterns and to determine the free energy of binding using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) approach. The extracted free energies from MM-PBSA calculations are in qualitative agreement with experimentally determined potencies. Furthermore, the interaction patterns seen in the receptor-GLP-1 complex simulations resemble previously reported binding interactions validating the simulations. Analysing the receptor-GLP-1 analogue complex simulations, we found that the major differences between the systems stem from FA interactions and positioning of acylation in the peptide. Hydrophobic interactions between the FA chain and a hydrophobic patch on the extracellular domain contribute significantly to the binding affinity. Acylation on Lys26 resulted in noticeably more interactions between the FA chain and the extracellular domain hydrophobic patch than found for acylation on Lys34 and Lys38, respectively. The presence of a charged linker between the peptide and FA chain can potentially stabilise the complex by forming hydrogen bonds to arginine residues in the linker region between the extracellular domain and the transmembrane domain. A molecular understanding of the fatty acid structure and its effect on binding provides important insights into designing acylated agonists for GLP-1R.Communicated by Ramaswamy H. Sarma.

Molecular insights into alginate ß-lactoglobulin A multivalencies-The foundation for their amorphous aggregates and coacervation.

For improved control of biomaterial property design, a better understanding of complex coacervation involving anionic polysaccharides and proteins is needed. Here, we address the initial steps in condensate formation of ß-lactoglobulin A (ß-LgA) with nine defined alginate oligosaccharides (AOSs) and describe their multivalent interactions in structural detail. Binding of AOSs containing four, five, or six uronic acid residues (UARs), either all mannuronate (M), all guluronate (G), or alternating M and G embodying the block structural components of alginates, was characterized by isothermal titration calorimetry, nuclear magnetic resonance spectroscopy (NMR), and molecular docking. ß-LgA was highly multivalent exhibiting binding stoichiometries decreasing from five to two AOSs with increasing degree of polymerization (DP) and similar affinities in the mid micromolar range. The different AOS binding sites on ß-LgA were identified by NMR chemical shift perturbation analyses and showed diverse compositions of charged, polar and hydrophobic residues. Distinct sites for the shorter AOSs merged to accommodate longer AOSs. The AOSs bound dynamically to ß-LgA, as concluded from saturation transfer difference and 1 H-ligand-targeted NMR analyses. Molecular docking using Glide within the Schrödinger suite 2016-1 revealed the orientation of AOSs to only vary slightly at the preferred ß-LgA binding site resulting in similar XP glide scores. The multivalency coupled with highly dynamic AOS binding with lack of confined conformations in the ß-LgA complexes may help explain the first steps toward disordered ß-LgA alginate coacervate structures.

Dynamics of Human Serum Transferrin in Varying Physicochemical Conditions Explored by Using Molecular Dynamics Simulations.

Conformational stability of human serum transferrin (Tf) at varying pH values and salt and excipient concentrations were investigated using molecular dynamics (MD) simulations, and the results are compared with previously published small-angle X-ray scattering (SAXS) experiments. SAXS study showed that at pH 5, Tf is predominantly present in a partially open (PO) form, and the factions of PO differ based on the physicochemical condition and drift toward the closed form (HO) as the pH increases. Tf is a bilobal glycoprotein that is composed of homologous halves termed the N- and C-lobes. The current study shows that the protonation of Y188 and K206 at pH 5 is the primary conformational drive into PO, which shifts toward the closed (HO) conformer as the pH increases. Furthermore, at pH 6.5, PO is unfavorable due to negative charge-charge repulsion at the N/C-lobe interface linker region causing increased hinge distance when compared to HO, which has favorable attractive electrostatic interactions in this region. Subsequently, the effect of salt concentration was studied at 70 and 140 mM NaCl. At 70 mM NaCl and pH 5, chloride ions bind strongly in the N-lobe iron-binding site, whereas these interactions are weak at pH 6.5. With increasing salt concentration at pH 5, the regions surrounding the N-lobe iron-binding site are saturated, and as a consequence, sodium and chloride ions accumulate into the bulk. Additionally, protein-excipient interactions were investigated. At pH 5, the excipients interact in similar loop regions, E89-T93, and D416-D420, located in the N- and C-lobes of the HO conformer, respectively. It is anticipated that interactions of additives in these two loop regions cause conformational changes that lead to iron-coordinating residues in the N-lobe to drift away from iron and thus drive HO to PO conversion. Furthermore, at pH 6.5 and 140 mM histidine, these interactions are negligible leading to the stabilization of HO.

pH- and concentration-dependent supramolecular assembly of a fungal defensin plectasin variant into helical non-amyloid fibrils.

Self-assembly and fibril formation play important roles in protein behaviour. Amyloid fibril formation is well-studied due to its role in neurodegenerative diseases and characterized by refolding of the protein into predominantly ß-sheet form. However, much less is known about the assembly of proteins into other types of supramolecular structures. Using cryo-electron microscopy at a resolution of 1.97 Å, we show that a triple-mutant of the anti-microbial peptide plectasin, PPI42, assembles into helical non-amyloid fibrils. The in vitro anti-microbial activity was determined and shown to be enhanced compared to the wildtype. Plectasin contains a cysteine-stabilised α-helix-ß-sheet structure, which remains intact upon fibril formation. Two protofilaments form a right-handed protein fibril. The fibril formation is reversible and follows sigmoidal kinetics with a pH- and concentration dependent equilibrium between soluble monomer and protein fibril. This high-resolution structure reveals that α/ß proteins can natively assemble into fibrils.

Investigation of the pH-dependent aggregation mechanisms of GCSF using low resolution protein characterization techniques and advanced molecular dynamics simulations.

Granulocyte-colony stimulating factor (GCSF) is a widely used therapeutic protein to treat neutropenia. GCSF has an increased propensity to aggregate if the pH is increased above 5.0. Although GCSF is very well experimentally characterized, the exact pH-dependent aggregation mechanism of GCSF is still under debate. This study aimed to model the complex pH-dependent aggregation behavior of GCSF using state-of-the-art simulation techniques. The conformational stability of GCSF was investigated by performing metadynamics simulations, while the protein-protein interactions were investigated using coarse-grained (CG) simulations of multiple GCSF monomers. The CG simulations were directly compared with small-angle X-ray (SAXS) data. The metadynamics simulations demonstrated that the orientations of Trp residues in GCSF are dependent on pH. The conformational change of Trp residues is due to the loss of Trp-His interactions at the physiological pH, which in turn may increase protein flexibility. The helical structure of GCSF was not affected by the pH conditions of the simulations. Our CG simulations indicate that at pH 4.0, the colloidal stability may be more important than the conformational stability of GCSF. The electrostatic potential surface and CG simulations suggested that the basic residues are mainly responsible for colloidal stability as deprotonation of these residues causes a reduction of the highly positively charged electrostatic barrier close to the aggregation-prone long loop regions.

Water Distribution on Protein Surface of the Lyophilized Proteins With Different Topography Studied by Molecular Dynamics Simulations.

Water play an important role in many structural and physicochemical properties of lyophilized proteins. Molecular dynamics simulations were employed to study the explicit water distributions on four structurally diversed proteins: insulin-like growth factor 1 (IGF1), immunoglobin G1 (IgG1), human serum albumin (HSA), and collagen. The MD simulations were combined with the literature data on water vapor sorption isotherms. To account for the heterogeneity of protein surface, the water molecules were classified into different groups according to the binding strengths. A mechanistic mathematical model was built to describe the type-II vapor sorption isotherms and successfully applied to all four model protein systems. Although commonly used Brunauer-Emmett-Teller (BET) theory has a good fitting to the experimental vapor sorption isotherms, the basic "monolayer" concept is not consistent with reality - covering too limited protein surface. Experimentally, several physicochemical properties did show a break point near the BET "monolayer" level. This study demonstrates that the water content threshold or BET "monolayer" is consistent with the onset of water cluster (n≥3) formation. Based on water distributions at different amino acid sidechains as well as the backbones, a simple formula was derived based on primary sequence and fractions of ordered secondary structures (i.e. alpha helix and beta sheet) to predict the BET "monolayer". We find that proteins with helical structural elements are more stable upon changes in water content compared to other protein architectures.

Self-Interactions of Two Monoclonal Antibodies: Small-Angle X-ray Scattering, Light Scattering, and Coarse-Grained Modeling.

Using light scattering (LS), small-angle X-ray scattering (SAXS), and coarse-grained Monte Carlo (MC) simulations, we studied the self-interactions of two monoclonal antibodies (mAbs), PPI03 and PPI13. With LS measurements, we obtained the osmotic second virial coefficient, B22, and the molecular weight, Mw, of the two mAbs, while with SAXS measurements, we studied the mAbs' self-interaction behavior in the high protein concentration regime up to 125 g/L. Through SAXS-derived coarse-grained representations of the mAbs, we performed MC simulations with either a one-protein or a two-protein model to predict B22. By comparing simulation and experimental results, we validated our models and obtained insights into the mAbs' self-interaction properties, highlighting the role of both ion binding and charged patches on the mAb surfaces. Our models provide useful information about mAbs' self-interaction properties and can assist the screening of conditions driving to colloidal stability.

Combination of high throughput and structural screening to assess protein stability - A screening perspective.

High throughput screening for measuring the stability of industrially relevant proteins and their variants is necessary for quality assessment in the development process. Advances in automation, measurement time and sample consumption for many techniques allow rapid measurements with minimal amount of protein. However, many methods include automated data analysis, potentially neglecting important aspects of the protein's behavior in certain conditions. In this study we implement small angle X-ray scattering (SAXS), typically not used to assess protein behavior in industrial screening, in a high throughput screening workflow to address problems of contradicting results and reproducibility among different high throughput methods. As a case study we use the lipases of Thermomyces lanuginosus and Rhizomucor miehei, widely used industrial biocatalysts. We show that even the initial analysis of the SAXS data without performing any time-consuming modelling provide valuable information on interparticle interactions. We conclude that recent advances in automation and data processing, have enabled SAXS to be used more widely as a tool to gain in-depth knowledge highly useful for protein formulation development. This is especially relevant in light of increasing accessibility to SAXS due to the commercial availability of benchtop instruments.

Electrostatics Drive Oligomerization and Aggregation of Human Interferon Alpha-2a.

Aggregation is a common phenomenon in the field of protein therapeutics and can lead to function loss or immunogenic patient responses. Two strategies are currently used to reduce aggregation: (1) finding a suitable formulation, which is labor-intensive and requires large protein quantities, or (2) engineering the protein, which requires extensive knowledge about the protein aggregation pathway. We present a biophysical characterization of the oligomerization and aggregation processes by Interferon alpha-2a (IFNα-2a), a protein drug with antiviral and immunomodulatory properties. This study combines experimental high throughput screening with detailed investigations by small-angle X-ray scattering and analytical ultracentrifugation. Metropolis Monte Carlo simulations are used to gain insight into the underlying intermolecular interactions. IFNα-2a forms soluble oligomers that are controlled by a fast pH and concentration-dependent equilibrium. Close to the isoelectric point of 6, IFNα-2a forms insoluble aggregates which can be prevented by adding salt. We show that monomer attraction is driven mainly by molecular anisotropic dipole-dipole interactions that increase with increasing pH. Repulsion is due to monopole-monopole interactions and depends on the charge of IFNα-2a. The study highlights how combining multiple methods helps to systematically dissect the molecular mechanisms driving oligomer formation and to design ultimately efficient strategies for preventing detrimental protein aggregation.

Physical constraints and functional plasticity of cellulases.

Enzyme reactions, both in Nature and technical applications, commonly occur at the interface of immiscible phases. Nevertheless, stringent descriptions of interfacial enzyme catalysis remain sparse, and this is partly due to a shortage of coherent experimental data to guide and assess such work. In this work, we produced and kinetically characterized 83 cellulases, which revealed a conspicuous linear free energy relationship (LFER) between the substrate binding strength and the activation barrier. The scaling occurred despite the investigated enzymes being structurally and mechanistically diverse. We suggest that the scaling reflects basic physical restrictions of the hydrolytic process and that evolutionary selection has condensed cellulase phenotypes near the line. One consequence of the LFER is that the activity of a cellulase can be estimated from its substrate binding strength, irrespectively of structural and mechanistic details, and this appears promising for in silico selection and design within this industrially important group of enzymes.

Binding Sites for Oligosaccharide Repeats from Lactic Acid Bacteria Exopolysaccharides on Bovine ß-Lactoglobulin Identified by NMR Spectroscopy.

Lactic acid bacterial exopolysaccharides (EPS) are used in the food industry to improve the stability and rheological properties of fermented dairy products. ß-Lactoglobulin (BLG), the dominant whey protein in bovine milk, is well known to bind small molecules such as fatty acids, vitamins, and flavors, and to interact with neutral and anionic polysaccharides used in food and pharmaceuticals. While sparse data are available on the affinity of EPS-milk protein interactions, structural information on BLG-EPS complexes, including the EPS binding sites, is completely lacking. Here, binding sites on BLG variant A (BLGA), for oligosaccharides prepared by mild acid hydrolysis of two EPS produced by Streptococcus thermophilus LY03 and Lactobacillus delbrueckii ssp. bulgaricus CNRZ 1187, respectively, are identified by NMR spectroscopy and supplemented by isothermal titration calorimetry (ITC) and molecular docking of complexes. Evidence of two binding sites (site 1 and site 2) on the surface of BLGA is achieved for both oligosaccharides (LY03-OS and 1187-OS) through NMR chemical shift perturbations, revealing multivalency of BLGA for EPS. The affinities of LY03-OS and 1187-OS for BLGA gave K D values in the mM range obtained by both NMR (pH 2.65) and ITC (pH 4.0). Molecular docking suggested that the BLGA and EPS complexes depend on hydrogen bonds and hydrophobic interactions. The findings provide insights into how BLGA engages structurally different EPS-derived oligosaccharides, which may facilitate the design of BLG-EPS complexation, of relevance for formulation of dairy products and improve understanding of BLGA coacervation.

Direct coordination of pterin to FeII enables neurotransmitter biosynthesis in the pterin-dependent hydroxylases.

The pterin-dependent nonheme iron enzymes hydroxylate aromatic amino acids to perform the biosynthesis of neurotransmitters to maintain proper brain function. These enzymes activate oxygen using a pterin cofactor and an aromatic amino acid substrate bound to the FeII active site to form a highly reactive FeIV = O species that initiates substrate oxidation. In this study, using tryptophan hydroxylase, we have kinetically generated a pre-FeIV = O intermediate and characterized its structure as a FeII-peroxy-pterin species using absorption, Mössbauer, resonance Raman, and nuclear resonance vibrational spectroscopies. From parallel characterization of the pterin cofactor and tryptophan substrate-bound ternary FeII active site before the O2 reaction (including magnetic circular dichroism spectroscopy), these studies both experimentally define the mechanism of FeIV = O formation and demonstrate that the carbonyl functional group on the pterin is directly coordinated to the FeII site in both the ternary complex and the peroxo intermediate. Reaction coordinate calculations predict a 14 kcal/mol reduction in the oxygen activation barrier due to the direct binding of the pterin carbonyl to the FeII site, as this interaction provides an orbital pathway for efficient electron transfer from the pterin cofactor to the iron center. This direct coordination of the pterin cofactor enables the biological function of the pterin-dependent hydroxylases and demonstrates a unified mechanism for oxygen activation by the cofactor-dependent nonheme iron enzymes.

Computing Cellulase Kinetics with a Two-Domain Linear Interaction Energy Approach.

While heterogeneous enzyme reactions play an essential role in both nature and green industries, computational predictions of their catalytic properties remain scarce. Recent experimental work demonstrated the applicability of the Sabatier principle for heterogeneous biocatalysis. This provides a simple relationship between binding strength and the catalytic rate and potentially opens a new way for inexpensive computational determination of kinetic parameters. However, broader implementation of this approach will require fast and reliable prediction of binding free energies of complex two-phase systems, and computational procedures for this are still elusive. Here, we propose a new framework for the assessment of the binding strengths of multidomain proteins, in general, and interfacial enzymes, in particular, based on an extended linear interaction energy (LIE) method. This two-domain LIE (2D-LIE) approach was successfully applied to predict binding and activation free energies of a diverse set of cellulases and resulted in robust models with high accuracy. Overall, our method provides a fast computational screening tool for cellulases that have not been experimentally characterized, and we posit that it may also be applicable to other heterogeneously acting biocatalysts.

Development of a fast screening method for selecting excipients in formulations using MD simulations, NMR and microscale thermophoresis.

Development of peptide therapeutics generally involves screening of excipients that inhibit peptide-peptide interactions, hence aggregation, and improve peptide stability. We used the therapeutic peptide plectasin to develop a fast screening method that combines microscale thermophoresis titration assays and molecular dynamics simulations to relatively rank the excipients with respect to binding affinity and to study key peptide-excipient interaction hotspots on a molecular level, respectively. Additionally, 1H-13C-HSQC NMR titration experiments were performed to validate the fast screening approach. The NMR results are in qualitative agreement with results from the fast screening method demonstrating that this approach can be reliably applied to other peptides and proteins as a fast screening method to relatively rank excipients and predict possible excipient binding sites.